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O ; ABSTRACT 

Context. Silicate minerals belong to the most abundant solids that form in cosmic environments. Their formation requires that a 
^ ■ sufficient number of oxygen atoms per silicon atom is freely available. For the standard cosmic element mixture this can usually be 

taken for granted, but this becomes a problem at the transition from the oxygen rich chemistry of M-stars to the carbon rich chemistry 
of C-stars. In the intermediate type S-stars most of the oxygen and carbon is consumed by formation of CO and SiO molecules, and 
left-over oxygen to build Si04-tetrahedrons in solids becomes scarce. Under such conditions SiO molecules from the gas phase may 
condense into solid SiO. The infrared absorption spectrum of solid SiO differs from that of normal silicates by the absence of Si-O-Si 
bending modes around 18 /jm while the absorption band due to Si-0 bond stretching modes at about 10 yum is present. Recently it has 
been reported (Hony et al. 2009) that exactly this particular characteristics is observed in a number of S-star spectra. 
Aims. We demonstrate that this observation may be explained by the formation of solid SiO as a major dust component at C/O 
00- abundance ratios close to unity. 

(— I _ Methods. The infrared absorption properties of solid SiO are determined by laboratory transmission measurements of thin films of 

^ ^ SiO produced by vapour deposition on a Si(l 1 1) wafer From the measured spectra the dielectric function of SiO is derived by using a 

I ■ Brendel-oscillator model, particularly suited for the representation of optical properties of amorphous materials. The results are used 

O ' in model calculations of radiative transfer in circumstellar dust shells with solid SiO dust in order to determine the spectral features 

^ I due to SiO dust. 

C/3 ■ Results. Comparison of synthetic and observed spectra shows that reasonable agreement is obtained between the main spectral 

■ characteristics of emission bands due to solid silicon monoxide and an emission band centred on lOyum, but without accompanying 

18/jm band, observed in some S-stars. We propose that solid SiO is the carrier material of this 10/jm spectral feature. 

Key words, circumstellar matter - stars: mass-loss - stars; chemically peculiar - stars: AGB and post-AGB 

>: 
o 

tH- ■ .. . 9.7 um and 18 um. The formation of silicon monoxide has been 

^. . Sihcon monoxide, SiO, is one of the most abundant and most ^^^^ however, to be an important intermediate step in the 

O : stable molecules encountered m space. It is found to be present formation of silicate dust for the following reasons: 
^-H , in many astronomical objects by observing its infrared molec- 

CN ■ ular bands. Usually (but not exclusively) it is found in objects _ The back-bone of the sihcate mineral structure are Si04- 

; where rather warm (T > 1000 K) gaseous material exists, for tetrahedra with four oxygen atoms attached to a silicon atom. 

^ . instance in circumstellar dust shells. At lower temperatures the The silicate minerals detected in space ai-e usually identified 

SiO molecules tend to associate with Mg, Fe, and O (available as by the stretching and bending vibrational modes of this par- 

H2O) from the gas-phase to form solid magnesium-iron-silicates ticular structure. 

^ . - the main components of cosmic dust - and some additional mi- _ On the molecular level there exists no species with compara- 

■ - ■ nor condensed phases bearing also Al and/or Ca. Alternatively ble structure as free gas-phase species, 
the SiO molecules may also condense to a solid material for their 

own, the silicon monoxide solid. Solid SiO is a well known ma- It is necessary, therefore, that the initial stages of the condensa- 

terial of significant technical importance, as it is widely used for tion of silicates involve some different kind of material. Because 

optical purposes as anti-reflection coating. Despite of its techni- of the high abundance of SiO molecules in matter with standai'd 

cal use, the properties of SiO are not well known, and its lattice cosmic elemental abundances it has been speculated since the 

structure remains enigmatic. early 1980s that condensation of SiO may be t he very initial st ep 

In astrophysics, soHd SiO has not yet been detected as a sep- of sil i cate form ation (Nuth III & Donn 1982j; iGail & SedlmavJ 

arate dust component. In the early discussions on the nature of Il986l Il998alltt iNuth III & Ferguson .2006i) . It is assumed that 

ISM dust it was taken into con sideration as possible ISM dust initially sihcon monoxide clusters are formed that serve as seed 

component (Dulev et al. 1978; see also Millar 1982', and ref- particles for growth of the silicates. 

erences therein), until amorphous silicates (olivine, pyroxene) There exists, however, the possibility that in a number of 

cases the silicon monoxide seeds finally grow to a dust species 

Send offprint requests to: gail@uni-heidelberg.de of their own. The formation of silicates requires the free avail- 
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ability of sufficient oxygen to build tlie Si04-tetrafiedrons in sili- 
cates. This is usually not a problem since the oxygen abundance 
in the standard cosmic element mixture is about twice the car- 
bon abun dance and about the sixteen-fold of the silicon abun- 
dance (cf. iL odders et al.ll2009l) . Even after formation of the ex- 
tremely stable CO and SiO molecules sufficient oxygen remains 
available to form the Si04-tetrahedra of silicate minerals from 
SiO. However, this does not hold during the whole lifetime of 
stars. During the evolution of low and intermediate mass stars 
on the asymptotic giant branch (AGB) the C/O abundance ra- 
tio stepwise increases once "third dredge-up" starts to operate. 
This drives in a number of steps the C/O abundance ratio to 
values exceeding unity, thereby passing through a stage where 
the C/O ratio is close to unity. During this stage the stars show 
the characteristic chemical peculiarities of the S-stars. Then al- 
most all oxygen is used up by the formation of CO and SiO, and 
with increasing C/O abundance ratio only a decreasing fraction 
of the SiO molecules can finally be converted into silicate dust 
because of the decreasing fraction of oxygen left over after SiO 
and CO formation. Since it is possible to form a solid material 
with chemical composition SiO by condensation of SiO vapour, 
there arises the possibility to form such kind of material in stellar 
outflows from S-stars where due to "third dredge-up" the avail- 
able oxygen becomes scarce. 

Recently some observational results for infrared emission 
from dust in S-stars have been published which show spectral 
characteristics that cannot be explained by the standard amor- 
phous silicate dust dHonv et al.ll2009t: ISacuto et al.ll2008 ). For 
these stars the emission feature at about 10 /im is as strong as 
usual, but the concomitant feature of silicate dust at 1 8 //m is 
very weak or even seem s to be missing at all. By inspection of 
the spectra published by IChen & Kwokl (1 1993b one finds some 
additional objects sharing this special property. In our opinion 
this property is a clear signature of the emission from solid sil- 
icon monoxide dust. In this material one has Si-O bonds as in 
silicates, which give rise to the stretching vibrations that are the 
origin of the broad lOyum feature of amorphous silicates, but one 
has no Si-O-Si bending modes as in Si04-tetrahedra which are 
the origin of the 18/im feature of silicates. Theoretically, this 
feature does not exist for solid silicon monoxide, but in practice, 
since solid silicon monoxide is prone to disintegrate into silicon 
nano-clusters and Si02 (e.g. iHohl et ani2003t Ivan Hapert et al.l 
l200l . the 18/im feature is not completely absent but weak. In 
S-stars with some oxygen left over after formation of CO and 
SiO molecules part of the SiO may form some silicate material 
while the bulk of the SiO molecules condense into solid silicon 
monoxide. This would give rise to an unusually weak 18//m fea- 
ture compared to the lOfim feature. 

In this paper we attempt to show that the observational find- 
ing of the existence of S-stars with strong IQfim emission fea- 
ture and weak or absent 18/im feature can be explained by con- 
densation of solid silicon monoxide in their outflows at C/O 
abundance ratios close to unity where chemically available oxy- 
gen to form silicates is lacking. For this purpose we calculate ra- 
diative transfer models of circumstellar dust shells with silicon 
monoxide dust and compare the resulting spectra with published 
spectra of S-stars. 

For these radiative transfer calculations one needs the dielec- 
tric function of solid silicon monoxide to calculate the extinction 
coefficient. Because this material is used for optical purposes the 
complex index of refraction is known from the middle infrared 
to the extreme UV and data are list ed, e.g., in Pal ik ( 198 5). More 
recent determinations are given bv lTazawa et al I ( |2006|) . During 
the course of our laboratory studies on materials of astrophys- 
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Fig. 1. Relative transmittance spectra of SiO on Si(lll) are 
shown for increasing film thicknesses up to 100 nm. The thick- 
ness change between two shown spectra is about 6.7 nm. The 
strong peak at 982 c m~' is attributed to the asymmetric Si-O 
stretching vibration dCachard et al.l [19711) . For such film thick- 
nesses, there is clearly no shift of that peak. Narrow features in 
the spectra are due to incomplete compensation of absorptions in 
the optical beam path of the purged spectrometer. SiO oscillator 
positions are indicated by dashed lines. 



ical interest (' Klevenzl 12009*) the far infrared optical properties 
have been determined with high accuracy. We briefly describe in 
this paper the measurements and their results. In particular the 
dielectric function is determined by fitting a Brendel-oscillator 
model (Brendel & Bormann 1992) to experimental results. This 
model is particularly suited to describe infrared optical prop- 
erties of amorphous materials in contrast to the widely used 
Lorentz-oscillator model that is better suited for crystalline ma- 
terials. The IQjjm feature of silicon monoxide is always diffuse 
like the corresponding feature of amorphous silicate materials. 
The lattice structure of solid silicon monoxide is not well known, 
but a crystalline phase seems not to exist. 

The plan of this paper is as follows. In Sect. |2] we briefly de- 
scribe our laboratory measurements on the optical properties of 
solid SiO and the results for the dielectric function. In Section |3] 
we give our arguments why solid SiO could form in outflows 
from S-stars. Section |4] describes the radiative transfer models 
and gives the results. 

2. Optical properties of SiO 

To study the infrared (IR) optical properties of SiO spectroscopic 
transmittance measurements were performed in situ on con- 
densed SiO films produced by evaporation under UHV condi- 
tions. The dielectric function of solid SiO is derived via Brendel- 
oscillator fits to the experimental IR spectra. 

2. 1 . Experimental setup 

The experiments were performed under UHV conditions (base 
pressure < 10"'*'mbar). IR spectra were taken in the range 
between 450 cm"' and 5000 cm"' with a Fourier-transform 
IR spectrometer (Bruker Tensor 27 with mercury-cadmium- 
telluride detector and N2 purged beam path). The IR spectral 
measurements were performed in transmittance geometry at nor- 
mal incidence of light. A floating zone Si(lll) wafer with di- 
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Fig. 2. Normal transmittance measurement of an evaporated SiO 
film (black solid line) with lOOnm thickness shown together 
with the best fit by a Brendel-model of the dielectric function 
as described in the text. Peak positions are indicated by dashed 
lines. 



mensions of 10 x lOmm^ and a thickness of 0.525 mm and high 
resistivity was used as transparent substrate which was cleaned 
by heating up to 1000°C in UHV to remove the natural oxide 
layer. Spectra were taken in situ, i.e., during the growth of the 
evaporated film on the substrate (spectral resolution of 4 cm '), 
and normalized to the spectrum of the bare Si(lll) taken directly 
before deposition. This method allows to monitor the whole 
growth process and to detect changes due to the film growth. 

Commercial SiO (Noah Technologies Corporation, Silicon 
Monoxide, SiO, 99.99% pure, -325 mesh, CAS no. 10097-28-6) 
was evaporated from a tantalum Knudsen cell heated by elec- 
tron bombardment. The background pressure measured during 
evaporation away from the SiO molecular beam did not exceed 
5 X 10"'° mbar. The SiO deposition rate (about 0.8 nm/min) was 
determined before and after the spectroscopic experiment with a 
quartz crystal microbalance with a relative error up to 5%. This 
error is caused by the uncertainty in the density of the evaporated 
film in the thickness calcul ation for which the den sity of bulk 
SiO (2.18 g/cm^) was taken (Hass & Salzberg' 1954. Fur ther ex- 
perimental details can be found in Klevenz et alT ( 2010ah . 



2.2. Results 

In Fig. [1] relative transmittance spectra of SiO with increasing 
film thicknesses are shown. The films were produced by SiO 
evaporation and subsequent condensation on a Si(lll) substrate 
hold at 300 K. Two broad features develop with growing film 
thickness. The stronger one at 982 cm"' is assigned to the asym- 
metric stretching mode of the oxygen atom in a Si-O -Si bridge 
against t he Si-neighbor atoms (Philipal97lt Chabal et al. ^200' 
[Oueenev et al.ll20()4l: ICachard et al.l(l97l ' 
[1984 1 



Lehmann et al.l l 19831 
The much weaker structure at approximately 713 cm"' 



corresponds to the Si stretching rnode, a mode with domina- 
ting Si displacement dLehmann et al.lll983i) . sometimes called 
"bending" mode in the literature. The third str ong p eak caused 
by the rocking mode at about 380 cm"' (Philipp(l97l]) was out of 
our measurement range. It is important to note that already from 
1 nm the IR spect ral feature s of the evaporated SiO films do not 
change anymore dKIevenz et al. 2010b). Accordingly, the IR di- 
electric function of the evaporated film already resembles that 
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Fig. 3. The dielectric function obtained from the best fit with four 
Brendel oscillators, as shown in Fig.|2] 



of bulk S iO (Hiortsber g & Granqvist|[T980l iTazawa et alj|2006l: 
iKlevenz et al.,.2010a.) . 



2.3. Dielectric function 

For thin dielectric films with the thickness d of the film much 
smaller than the wavelength A {d <^ A), approximate formu- 
las for the transmittance ca n be derived based on the Fresnel 
equations (iBerrem anl ll963h . For a thin layer on a thick (non- 
int erfering) substra t e the normal transmittance {(p - 0°) is given 
by (lLehmannlll988l:lTeschner & Hubneilll990h 



1 - 



2d^ 



1 + 



Im (sf) , 



(1) 



with 0) - Ittc/A, £s as the dielectric function of the substrate and 
fif as the dielectric function of the film. In transmittance mea- 
surements under normal incidence {<p - 0°) therefore only the 
transverse optical (TO) modes with frequency at the maximum 
of the imaginary part of the dielectric function s and with a dy- 
namic dipole moment parallel to the surface can be observed. 

To determine the dielectric function of SiO in the IR from 
such thin film measurements an appropriate model for the dielec- 
tric function has to be used. A simple Lorentz-oscillator model 
can not be applied in our case due t o the am orphous character of 
the film. The model introduced bv iBrende l & Bormann (1992) 
accounts for the amorphous structure of a material by assum- 
ing that the different IR modes can be represented by Lorentz- 
oscillators, but with randomly shifted resonance frequencies that 
are distributed according to a Gaussian probability distribution. 
The dielectric function in the IR therefore is assumed in this 
model to be given by the relation 



g(<jj) 



IR 



dze 



z — ci)-'- + iyjto 



(2) 



It consists of a dielectric background Soo and distributions 
of Lorentz oscillators with resonance frequencies wo y, damping 
constants -yj, plasma frequencies cuipj, and with standard devia- 
tions (Tj of the Gaussian probability distributions. 

There is a close relationship of the Brendel oscillator 
model to the well known Voigt function for the profiles of 
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Table 1. Constants of the Brendel oscillator model for the di- 
electric function of amorphous materials, Eq. (|2|l, for our best fit 
to our experimental results on sohd SiO. All quantities are given 
in cm"' units. 



Osc. 


CDq 


CUp 


tr 


r 


1 


1101 


305 


44.9 


4 


2 


982 


699 


57.4 


4 


3 


713 


298 


73.4 


4 


4 


380 


461 


127.1 


4 




3.61 









velocity broadened spectral lines. This is briefly discussed in 
AppendixiBl 

The last spectrum of Fig. [T]was fitted with such a Brendel 
dielectric func tion using four oscillators, similar to a proce- 
dure for SiO-, (lKirklll9"88h . which results in very good agree- 
ment (see Fig. |2]l. We performed the spectral fits in the range 
between 500 cm' ' ' and 1500 cm"' using the software package 
SCOUT dTheissI [2009). For the high frequency limit of 
the dielectric function a value of 3.61 was taken from litera- 
ture (iHiortsberg & GranqvisjlT980h . Since the oscillator damp- 
ing constant y cannot properly be dete rmined from fits to vi bra- 
tion spectra of disordered solids (cf . Ilshikawa et alj l2000l) we 
fixed the value of y to the resolution of our measurement of 
4 cm"' . Nearl y the same kind of procedure is already reported in 
the hterature (Ilshikawa et alj|2000t iBrendel & BormannI [19921 



lGrosse etal.]ll986l:lNaimanetaljri985h . The result for the di- 
electric function of SiO in the IR is shown in Fig. [3] and our best 
fit for the oscillator parameters is given in Table [1] 

These results are obtained for a substrate temperature of 
300 K. Other substrate temperatures than 300 K result in slightly 
shifted peak positions. At a substrate temperature of 93 K the po- 
sition of the asymmetric stretching mode is at 969 cm"', while 
for 300 K it is at 982 cm"'. For tempe ratures hi gher than 300 K 
the peak shifts to higher wavenumbers (lKlevenzi200 9). The shift 
to higher wavenumbers can be explained by a decomposition 
into Si and SiOa that accelerates with temperature. The effects 
of the substrate temperature on the condensed SiO film and the 
dielectric function of SiO will be discussed in detail in a forth- 
coming paper. In our present communication the dielectric func- 
tion obtained at 300 K is used. 

3. Astrophysical applications 

In this section we discuss the possibility that solid silicon 
monoxide might be an important dust species in the dust shells 
of S-stars. 



3.1. S-stars 

Low and intermediate mass stars from the range of initial masses 
below X 8 M0 evolve till the very end of their lives to the ther- 
mally pulsing Asymptotic Giant Branch (TP-AGB) where they 
are composed of an electron-degenerated core of carbon and 
some oxygen, resulting from He-burning, an overlying layer of 
mainly He, resulting from H-burning via the CNO-cycle, and an 
enormously extended hydrogen rich envelope. On the TP-AGB 
the stars alternatively burn either H in a shell source at the inter- 
face between the H-rich envelope and the He layer for a period 
of several thousand years, or He in a shell-source at the interface 
between the He layer and the carbon-oxygen core for a period 



of about 200 yr. During the short burning phase of He (the ther- 
mal pulse) the He layer is for part of the time fully convective 
from bottom to top and distributes freshly synthesized carbon 
from the He-burning shell at the bottom of the He layer over the 
whole layer Somewhat later in the pulse phase, the convection 
zone of the fully convective hydrogen envelope briefly dips into 
the upper parts of the He layer and mixes some carbon-rich mate- 
rial into the outer envelope. As a result, after each thermal pulse 
the carbon abundance of the envelope increases stepwise. The 
oxygen abundance in the envelope remains almost unchanged 
by this process because only very little oxygen is synthesized by 
He-burning in TP-AGB stars. Thus, the carbon to oxygen abun- 
dance ratio C/0 in the envelope changes from its initial value of 
~ 0.5, corresponding to t he C/O ratio of the c osmic standard el- 
ement mixture (see, e.g.. iLodders et al]|2009l) . to C/O ratios far 
exceeding a value of unity. 

For most of the stars on the TP-AGB the C/O abundance 
ratio falls — after mixing carbon from the core into the en- 
velope during a thermal pulse — for one or two inter-pulse 
phases into the critical range between about 0.9 and 1 .0 where 
the chemical composition of the material in the stellar photo- 
sphere dramatically changes from an oxygen-compound domi- 
nated composition for a C/O abundance ratio below 0.9 and a 
hydrocarbon-dominated composition for a C/O abundance ratio 
exceeding 1.0. This switching of the chemical composition is 
brought about by the extraordinary stability of the CO molecule 
that almost completely consumes by its formation the less abun- 
dant of the two elements C and O. Only part of the TP-AGB 
stars leap over the critical range of C/O abundance ratios during 
a single thermal pulse. 

In the intermediate abundance range of C/O ratios between 
0.9 and 1 .0 the stellar spectrum is dominated by molecular bands 
of some low-abundance elements that are neither seen in M-Stars 
(C/O ratio < 0.9) nor in C-stars (C/O ratio > 1 .0). These stars are 
the so called S-stars. They are much less abundant than M-stars 
or C-stars since stars on the TP-AGB suff'er numerous thermal 
pulses, but only during one or two inter-pulse phases they are 
S-stars. 

For these S-stars the material that outflows from the stellar 
surface neither can form the silicate dust that is seen in circum- 
stellar dust shells around M-stars, nor the graphite and SiC dust 
seen in dust shells around C-stars. This is because the excess 
oxygen over carbon, that is not bound in CO molecules, is too 
rare to combine with all of the silicon to the Si04-tetrahedrons 
that form the back-bone of silicate minerals. At the same time 
no excess carbon is available to form solid carbon. 

The details of the chemi stry of dust formation in outflows 
from S-stars is discussed in Ferrarotti & Gail' (■2002). Because of 
the very high bond energy of the SiO molecule — not as high 
as for CO, but also exceptionally high — the silicon is bound in 
this molecule and the excess of the oxygen not bound in SiO 
and CO forms H2O. Some quantities of silicate dust may be 
formed in outflows from S-stars if the C/O abundance ratio is 
not too close to unity, because then SiO molecules react with the 
left over H2O and with Mg (and/or Fe). Indeed, weak emission 
bands showing the characteristic two peaks around 9.7 fim and 
1 8 jum are seen in the infrared spectrum of a numbe r of S-stars 
(iChen & Kwok 1993; Llovd Eva ns & Little-Mareninl'l999h . It is 
argued by Ferrarotti & Gail (2002^ that also iron dust may be 
formed and may be an abundant dust species for S-stars. This 
dust would be hard to detect observationally for optically thin 
dust shells, since its smooth and featureless extinction would act 
as a more-or-less gray opacity source. For optically thick shells 
it is likely that this is mistaken to be carbon dust because of al- 
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Fig. 4. Stability limits against vaporisation of some refractory 
compounds formed from an oxygen-rich element mixture with 
cosmic standard element abundances, and of solid SiO. Dashed 
line: typical pressure-temperature combinations in an outflow- 
ing gas element for a stellar wind model with stationary outflow 
and M = 5 X lO"''M0yr"' (here for an M-star). Most part of 
the dust condensation occurs around the knee where the velocity 
turns from sub-sonic (horizontal part of the curve) to highly su- 
personic (steep descending part of the curve) outflow. Pressures 
are typically 10 bar in the region of main dust growth. 



most indistinguishable infrared extinction properties of iron and 
carbon dust. 



3.2. I/Vhy so//d silicon monoxide? 

If the abundance of excess oxygen not bound in CO and SiO 
tends to zero for C/O abundance ratios approaching unity, there 
exists the possibility that SiO molecules start to condense as a 
solid of their own. This has not yet seriously be taken into ac- 
count as a dust component in circumstellar dust shells, though it 
was discussed several times whether it could be the first conden- 
sate to form the necessary seed particles for silicate dust forma- 
tion (Nuth III & Donn 1982t lOail & Sedlmavrl 119861 il998allbl: 
iNuth III & Fergusoa2 006). This is a crucial question for cosmic 
dust formation since it is not possible to form the necessary seed 
particles for silicate dust growth from a material with the compo- 
sition and lattice structure of one of the many silicate minerals. 

The main reason why solid SiO was not considered as one of 
the major dust materials in circumstellar shells was the seem- 
ingly high vapour pressure from which one predicts a rather 
low condensation temperature under the low pressure conditions 
of circumstellar dust shells. New det ermination s of the vapour 
pressure of solid SiO (iFerguson & Nut h III 2008t|Ki even zl2Q09h 
have now shown that the older measurement were seriously 
in error and s trongly overestimated the vapour pressure. The 
new results of Ferguson & Nuth III ( 2008) and our own results 
dUevenz 2009) can be approximated by a vapour pressure for- 
mula for the equilibrium pressure of SiO molecules over solid 
silicon monoxide 



/ 42648 
p = exp H- 33.848 



(3) 



where the pressure is in units of dyncm ^ (as usual in astro- 
physics). Our results dKleven 1^009) extend the range of mea- 
sured vapour pressures down to the pressure range relevant for 



Table 2. Some data used for calculating the opacity of the dust 
species considered in the model calculations and the correspond- 
ing key elements 



species 


Ad 


Pd 


El. 


£ 




SiO 


44.09 


2.13 


Si 


3.55 10- 


-5 


iron 


55.85 


7.87 


Fe 


3.16 10- 


-5 


olivine 


172.2 


3.81 


Si 


3.55 10- 


-5 


pyroxene 


116.2 


3.61 


Si 


3.55 10- 


-5 



circumstellar dust shells. No extrapolation of vapour pressures 
measured at higher temperatures in the laboratory to signifi- 
cantly lower temperatures is required. The details will be dis- 
cussed in a separate paper 

Figure |4] shows the resulting stability limit against evapora- 
tion in chemical equilibrium for the astrophysically important 
silicates, for solid iron, and for solid SiO, calculated by using 
the new results for th e vapour pr essure of solid SiO and ther- 
mochemical data from lBarinI (1199 5) for the othe r materials, for 
the case of standard cosmi c element abundances dAsplund et al.l 
120051: iLodders et al]|2009l) . At C/O abundance ratios close to 
unity no silicates can be formed in S-stars because of a lack of 
freely available oxygen to convert a significant fraction of the 
SiO molecules from the gas phase into the Si04-tetrahedrons in 
silicate minerals (e.g., for a C/O ratio of 0.95 only about 20% 
of the Si can condense as silicate dust). Then the SiO molecules 
may condense instead into solid SiO. At a typical pressure of 
p = 10""^ bar in the dust condensation layer of circumstellar 
dust shells (see Fig. HJ the SiO would become stable against 
evaporation at about 800 K. This is high enough that such dust 
could be observed as warm dust in a circumstellar dust shell. 
Hence solid SiO is a candidate dust material to be observed in 
S-starsQ 

Whether this material really forms can be tested observation- 
ally. The absorption properties of solid SiO are sufficiently dif- 
ferent from normal silicate materials that they can be discrimi- 
nated by their spectral features. Amorphous olivine and pyrox- 
ene both have two strong absorption features centred on 9.7 //m 
and 18/im, that are usually observed as emission features in 
the infrared spectra from circumstellar dust shells. For SiO the 
strong feature at 18//m is missing because this corresponds to 
rocking modes in the Si04-tetrahedron, but a strong feature at 
about IQfim is present because this originates from Si-O bond 
stretching vibrations. Hence, the presence of a strong lOjt/m fea- 
ture without a strong 18yum feature is diagnostic for solid SiO. 

These characterist ics have, indeed, been observed for a num- 
ber of S-star spectra dHonv et al.ll2009 ) and one may speculate 
that they result from solid SiO as a major dust component in 
such dust shells. In order to test this hypothesis we have per- 
formed radiative transfer calculations using our new results for 
the dielectric function of SiO. 



3.3. Growtti of silicon monoxide grains 

Now we turn to the question whether growth of silicon monoxide 
is kinetically possible. Consider the simple approximation that a 



' A lack of oxygen to form normal silicates may also be encountered 
in massive supergiants if the outer layers are peeled off by mass-loss 
and material that has burned H via the CNO-cycle appears at the sur- 
face. Then the O abundance may drop to values of the order of the Si 
abundance and solid SiO may form also in the outflows from such stars. 
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spherical dust grain co-moving with the outflowing gas collects 
SiO molecules from the gas phase. For simplicity we neglect 
evaporation. The equation of growth of radius a for a single grain 
is 



da 

Vexp— = VdaWgr Vgr<P(f/drift) 



(4) 



v'expdf the radial coordinate instead of 



if we introduce by dr 

time as independent variable, Vexp being the outflow velocity. 

The meaning of the quantities in the growth equation is as 
follows: Vd is the volume of one chemical formula unit in the 
condensed phase 



Pd 



(5) 



Here Ad is the atomic weight of the condensed phase and pd its 
mass-density. Numerical values are given in Table|2] 

The quantity a is the growth coefficient. This seems to be 
first measured by Gunther ( 1958) who found a rather small value 
of g 4 X 10'^ in th e temperature range 1 200 K . . . 1 500 K. 
iRocabois et £0(119921) made a new determination of the vapour 
pressure and of the evaporation coefficient a of SiO and found 
that their results for a can be approximated by 



a(r) = 0.1687 ■ 



2.909 X lO"'* r + 1 .373 x 10"^ 



(6) 



in the temperature range 1 175 K to 1410K. Essenti a lly th e 
same results for a are found by iFerguson & Nuth nil (I2008h . 
Extrapolating these results down to a temperature of 800 K re- 
sults in Of a: 0.025. This value will be used in the following con- 
siderations. 

The quantity «gr is the particle density of the growth species. 
If we concentrate on the initial growth phase where not yet much 
of the growth species is consumed from the gas phase, we can 
approximate «gr by 



M 



Wgr — 



1 .4/77 H 4;rr^v, 



(7) 



exp 



where tgr is the element abundance of the element (Si in our 
case) that determines the gas-phase abundance of the growth 
species. The second quantity on the r.h.s. is the density of H 
nuchdes in a spherically symmetric, stationary outflow with ve- 
locity Vexp and mass-loss rate M. 
The quantity 



kT 



(8) 



is the thermal velocity of the growth species (SiO molecules in 
our case) and A„i- their molecular weight, while 



Cl)(f/): 



1 + 



16vl 



(9) 



is a correction factor for particle drift relative to the gas with 
relative velocity U. If U » 4vgi the particle drift is supersonic 
with respect to the thermal velocity of the growth species (not 
the carrier gas!); then Vgr<i>{U) ~ jU. 

Typical values of all these quantities required for calculating 
SiO growth are given in Table[3] 

For the following we assume the outflow velocity Vexp in the 
region where dust grows to be constant. This is not completely 
correct since dust condensation results in an acceleration of the 



Table 3. Data used for calculation of silicon monoxide growth 



Quantity 


unit 


value 


Vd 


cm^ 


3.46 X 10"^^ 


A r 




44.09 






3.55 X 10-^ 


Ed 




5 X 10-" 


a 




0.025 


^max 


yum 


0.0837 


Tc 


K 


800 


Vth 


cms"' 


1.54 X lO"* 




cm 


3.2 X 10'3 



wind, but the velocity increase is moderate if we consider S- 
stars. It increases from a velocity of about the sonic velocity 
(2kms-' - 3kms-') at the inner edge of the dust shell to the 
final out flow velocity that is at most lOkms"' for S-stars (e.g. 
Ramstedt etal ] l2006l) . Then integrating from the inner radius 
of the dust shell, where dust commences to grow, to infinity we 
obtain 



, - uq + Via Vgi-<1) fgi- 



M 



(10) 



■^exp 



where oq is the radius of the seed nuclei for dust growth and 
Qoo the grain radius at infinity. This equation requires that the 
grain radius Aoo remains smaller than the maximum radius amax, 
attained if all condensable material is condensed, because the 
consumption of the growth species is neglected in Eq. (|4|. The 
maximum possible radius to which a particle may grow in the 
outflow is given by 



47r 3 



(11) 



if ed is the number of dust grains per hydrogen nucleus. This 
quantity is not precisely known; we assume a value of 5 x lO"'"' 
that is a typical value found by observations (iKnappil985i) . 
We can now write 



^oo — ^0 ^max 



M 



with 



Vd a Vgr<l>% 



(12) 



(13) 



The mass-loss rate M has to be smaller than M^r in order that our 
assumption is valid, that there is no strong depletion of the gas 
phase from condensable material. 

For applying this to silicon monoxide condensation we have 
to specify the values of Vexp, Rc and All other quantities are 
given in Table |3] For the outflow velocity we assume a value 
of Vexp - 3kms-', slightly higher than the sound velocity of 
the gas, since most of the grain growth proceeds before the gas 
is accelerated to highly supersonic outflow velocities; otherwise 
further growth is suppressed by rapid dilution of the gas. For R^ 
we assume a value /?c = 3 x lO''' cm, corresponding to about 
8 stellar radii, see Table |4] For the drift velocity we assume the 
same value as for Vgxp since typical drift velocities usually are of 
the order of the sonic velocity of the gas, except for very high 
mass-loss rates where they are much less. With these estimated 
values we obtain 



= 4.5x 10"^ Moyr- 



(14) 
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This shows that significant condensation of sohd silicon monox- 
ide is kinetically possible, at least for mass-loss rates exceeding 
about 10"'' M0. 

The result for M^r.sii depends on the rather uncertain value of 
a that is determined from extrapolating laboratory measured val- 
ues (which decrease with decreasing temperature) to much lower 
temperature. The value of 0.025 at 800 K determined this way 
may be too low. In this case Mcr,sii would be overestimated and 
growth of solid silicon monoxide may be even more favourable 
as in the estimate above. 

4. Radiative transfer model for dust shells 

4.1. Opacity 

The spectral energy distribution of the emission from circum- 
stellar dust shells is determined by the composition of the dust 
mixture and the properties of the different dust species. Here we 
concentrate on the special question whether the unusual feature 
at 10/im detected in a number of S-stars could result from solid 
SiO. Since this feature is the most prominen t dust feature in th e 
objects where it has been detected so far jHonv et al.l l2009h . 
it appears also to be the most abundant dust species in these 
objects, apart from possible contributions to opacity from dust 
species with structure-less pure continuous extinction like, e.g., 
iron dust grains. Therefore we concentrate on models with solid 
silicon monoxide as the sole dust component. 

In our model calculations we assume the opacity in the shell 
to be completely determined by the dust components; no con- 
tribution of the gas phase is considered. The dust particles are 
assumed to be small spheres of radius a. The absorption and 
scattering coefficients of the different dust components, charac- 
terized by an index /, then are given by 



,,abs,sc 



4 (1 H-4eHe)Pd,/ 



(15) 



The quantity A, is the molecular weight corresponding to the 
chemical formula of the condensed phase, e,- is the elemental 
abundance of a characteristic key element that is required to 
form the condensed phase, and is the fraction of the atoms 
of this element that is really condensed into the solid phase, p^j 
is the mass-density of the solid. The quantities (^^^•'''^ are the 
absorption and scattering efficiencies, divided by the particle ra- 
dius a. These qua ntities are calculated by means of Mie-theory 
(see lBohren & Huff man 1983) from the dielectric function of the 
dust material. The basic data used in the calculation of dust ex- 
tinction are given in Table |2] 

With respect to the particle shape we consider only spheri- 
cal grains. If the grains would be non-spherical, the absorption 
band profiles of dust grains would be distorted and their cen- 
tre would be shifted compared to the case of spherical grains. 
However, SiO grains formed in laboratory condensation experi- 
ments from the gas ph ase generally seem to be almost spherical 
dKamitsuii et alJl2004 their Fig. la). Since there is no obvious 
reason why this should be different for the tiny grains condensed 
in a circumstellar shell and for nanometre sized laboratory con- 
densed grains, this leaves no room for introducing a distribution 
of grain shapes like, e.g., the popular continuo us distribution of 
ellipsoids (see, e.g.. iBohren & Huffmanl[T983h . in order to ac- 
count for possible non-sphericity effects. 

With respect to the particle size we assume the canon- 
ical value of a = 0.1 yum that se ems to be typical for 
dust grains in circumstella r shells (e.g. iKruszewski et al.lll968t 
ISerkowski & Shawlll200Ih . Without calculation of a complete 
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Fig. 5. Absorption and scattering coefficient per unit mass of 
some dust species in a dusty gas with standard cosmic ele- 
ment abundances. Full line: solid SiO. Dashed line: amorphous 
olivine: Dotted line: amorphous pyroxene. 



model for the dust shell including hydrodynamics and dust 
growth the grain radius cannot be specified more precisely. 
Fortunately, for the far infrared spectral region the opacity does 
not depend on the grain size a since circumstellar grains are in 
any case much smaller than the wavelength A in the far infrared. 
Only for the absorption of stellar radiation there may be a slight 
dependency on grain size since in the optical and near UV spec- 
tral region the grains may satisfy only marginally the condition 
Ina/A <s: 1 that particles can be considered as small. 

For the purpose of model calculations of dust shells of cool 
giant stars one needs optical data for the dust species at least 
in the wavelength range from about 0.4yum to about 100 //m in 
order to cover (i) the wavelength range of the stellar radiation 
field from the ultraviolet to the mid infrared, and (ii) to cover 
the wavelength range of dust emission from the shell that ranges 
from the near infrared to the sub-mm region. Our own measure- 
ments discussed in Sect. |2] cover only the mid to far infrared 
wavelength range. Therefore we have to augment our data by 
data from other sources to cover also the spectral region from 
near infrared to UV. 

For the visible to far-ultraviolet s pectra l regions {A < 0.8 yum) 
we use data for solid SiO listed in iPalikI (11985 ). The real part 
of the complex index of refraction of our infrared extinction 
data and the data given in Palik nicely fit in the near infrared 
where they overlap. The imaginary part becomes very small in 
the wavelength region below 8 fim in our results and it is omit- 
ted in the range between 8 and 0.8 /im in the listing of Palik. In 
this region the absorption coefficient of solid SiO would become 
many orders of magnitude smaller than the scattering coefficient 
(or would vanish at all if data of Palik are used). In order to avoid 
numerical problems in solving the radiative transfer problem we 
arbitrarily increase the absorption in the wavelength range 0.8 to 
8 ijm by assuming that the dust material contains tiny inclusions 
of pure iron particulates with a volume filling factor of fy - 10"^ 
and calculating the average dielectric function (e) for this com- 
posite material from the M axwell-Garnett mixing rule (cf., e.g., 
IBohren & Hu"ffmai3[T98l 



(1 -/v)em+/j6chi 
l-/v+/vy8 



(16) 
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Table 4. Basic parameters of the dust-shell models 



Star 

Wind 
Dust 
Shell 



^exp 

M 
^abs 

/sio 
R. 



2 700 
1 X 10* 
3.18 X 10" 
= 534 

10 

1 X 10^^ — 3 X 10-5 

amorphous SiO 
0.5 

depends on M 
1 X 10' 
750 



K 

Lo 
cm 

Ro 

km s"' 
Moa-' 



K 



where 



3 err, 



+ 2em 



(17) 



Here £^ is the bulk dielectric function of the matrix material 
(solid SiO in our case) and fine the bulk dielectric function of 
the material forming the inclusions (solid iron in our case). 

This kind of additional absorption by impurities has no influ- 
ence on the calculated temperature of dust grains since this is de- 
termined by the strong absorption at A < 0.8 fim and also has no 
influence on the calculated radiation field. It is unlikely that any 
real dust material is pure and ideally transparent, and nano-sized 
iron inclusions are a fundamental constituents of, e.g., GEMS 
(e.g.lBradley 2003). 

The resulting mass-absorption and scattering coefficients 
of small particles of solid SiO are shown in Fig. |5] They 
are calculated according to Eq. ( fTSl l using the data given in 
Table |2l and calculating optical constants for SiO according to 
the Brendel-model for the dielectric function of amorphous ma- 
terials, Eq. (|2]|, using constants from Table [1] and following the 
pro cedure described above. 

iBrendel & Borman nl (fT992h have shown that the integral in 
Eq. (|2]i can be solved in terms of the complex probability func- 
tions and Gaussians. We did not make use of this in our opacity 
calculations but preferred to evaluate the integrals numerically 
since evaluation of the analytic expressions seems not to offer 
any particular advantage due to their complicated nature. 

For comparison also mass-absorption and scattering coeffi- 
cients of amorphous olivine and pyroxene are shownQ that are 
the main dust species in circumstellar dust of M-Stars. It is seen 
that solid SiO indeed shows a prominent absorption feature cen- 
tred around about 10 //m and only some weak structure on the 
long wavelength side. 

4.2. Model calculations 

Using the above opacity description we have calculated radia- 
tive transfer models of circumstellar dust shells with solid SiO 
as dust component for a range of mass-loss rates in order to de- 
termine the spectral features that would result from solid sil- 
icon monoxide. The essentials of the model of the dust shell 
and of the radiative transfer calculations are briefly described 
in Appendix lAl 

A fundamental parameter of the model is the location of the 
inner boundary of the dust shell, that depends on the mech- 



^ Data for silicates with x = 0.7, from the Jena-St. Petersburg data 
basis, accessible via: http:/www. mpia-hd.mpg.de/HJPDOC/ 



anism by which the dust forms from the gas-phase by nucle- 
ation and subsequent grain growth. About the details of these 
processes almost nothing is known and we take recourse to the 
simple assumption that the dust suddenly appears at some pre- 
scribed dust temperature and the fraction /sio of material con- 
densed into the solid phase is constant across the dust shell. This 
is the kind of approximation on which most existing models of 
circumstellar dust shells are based on. For the temperature Tc we 
assume for silicon monoxide dust a value of 750 K because the 
temperature where solid SiO becomes stable against vaporisa- 
tion under conditions in circumstellar dust shells is about 800 K, 
see Sect. 13.21 and because in dust shells this temperature proba- 
bly is somewhat lower if some super-cooling of the outflowing 
material should be required f or the onset of nuc l eation and con- 
densation. The old results of iGail & Sedlmav^ (1198 6') seem to 
indicate a temperature for onset of SiO condensation as low as 
600 K, but this cannot be maintained since that calculation is 
based on older vapour pressure measurements of solid SiO for 
which it is n ow known that they h ave su bstantially to be revised 
downwards (i Ferguson & Nuth Iljllooi lKlevenzll2009h . which 
increases the condensation temperature. 

For the basic parameters T^ff, L*, Vgxp we use representa- 
tive average values, since we do not intend to model individual 
stars but to perform only some explorative calculations. For the 
effective temperature of the central star we chose a value of 
2 700 K. This seems t o be a representative value for S stars (e.g. 
iKerschbaum & Hronl [T998: Kerschbaum 1999). The luminosity 
is chosen to be L = 10"* L0, which is typical for stars at the up- 
permost part of the TP-AGB, where one expects AGB-stars to 
pass through the stage of S-stars. Observed values given in the 
literature for stars with substantial mass-loss rate are of this or- 
der of magnitude or slightly less (e.g. .Joris sen & Knapp 19981 
iGroenewe^en & de Jonglll998t iRamstedt et al.i 2006).'The ob- 
served expansion velocities Vexp of S-stars vary between a few 
kms ' and abo ut 20kms~' and are typically between 5kms~' 
and 10 kms ' (*Groen ewegen & de Jong! Il998l : IRamstedt et al.l 
2006). We use a value of lOkms"'. The observed mass-loss 
rates of S-stars vary over a large range between values as low 
as a few times I0~^ Mc?, yr~' to s everal times 10~^ M(r, yr~' (e.g. 
Gro enewegen & de Jon g 1998; Ramstedt et al■ll2006^ ■ We con- 
sider M as free parameter and calculate models for different val- 
ues of the mass-loss rate. 

The basic model parameters of the calculations are shown in 
Table m 



4.3. Results for spectral energy distribution 

Models are calculated for mass-loss rates ranging from M = 
1 X lO"''M0yr"' to M = 3 x lO^^Moyr"' with solid silicon 
monoxide as the sole opacity source. Figure |6] shows the result- 
ing spectral energy distribution of the radiation emitted by the 
star and its dust shell. Since we assumed for the star a black- 
body radiation field, the strong structure of cool stellar spec- 
tra resulting from molecular absorption bands is missing in our 
models. This simplifies the identification of the dust features in 
the spectrum. Emission from warm dust is seen only at wave- 
lengths longer than 8 fim since the absorption coefficient of solid 
SiO is small at shorter wavelengths, see Fig. |5] The prominent 
absorption band around lOyum is clearly seen as a strong emis- 
sion feature in all of the spectra. Some structure is seen in the 
spectrum at longer wavelengths but there is no strong emission 
feature around 1 8 fim. 
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Fig. 6. Synthetic spectra of models of circumstellar dust shells 
with solid silicon monoxide dust particles for different mass-loss 
rates from M = 1 x 10"'' M© yr"' to M = 3 x 10"^ yr"'. 
Dotted line: stellar radiaton. 



This is similar to what is found in a number of S-star spectra. 
Figure [T] shows a more close-up view of our synthetic spectra 
in the wavelength region between 5 yum and 25 jim. For clarity 
the black body radiation field of the star is subtracted from the 
synthetic spectra and the resulting emission profiles are shifted 
relative to each other These emission profiles reflect the peculiar 
characteristics of the emission coefficient of amorphous solid sil- 
icon monoxide: A strong emission band centred around lOyum 
and only weak to negligible structure in the 15/im to 25 
wavelength region. Figure Q also shows one of the two different 
profiles of the emission bands around lOyum observed in S-stars 




M = 1- 10 



25 



Fig. 7. Emitted energy flux from star and dust shell, subtracted 
by black-body emission from the star (lines with crosses) for dif- 
ferent mass-loss rates (i n units M(T| yr"'), a nd averaged emission 
profile of S-stars from Honv et al.l (l2009l) . The S-star line pro- 
file is scaled such that the maxima of observed and calculated 
profiles coincide; for clarity, the spectra are shifted by arbitrary 
amounts. 



where no or only a very wea k 18/im feature is observed. The 
data are taken from Fig. 6 of iHonv et al.l (1200 9^. The profiles 
in Fig . Ill corresponds to the lower one in Fig. 6 of Hony et all 
(l2009h . This profile is scaled such that the peak value coincides 
with that of the synthetic emission pr ofiles in our models. 

The upper profile of Fig. 6 of Ho nv et all (|2009) is not con- 
sidered since it peaks at a significantly higher wavelength. It may 
result from a different kind of material or a complex superposi- 
tion of emission bands from a number of materials. 



4.4. Discussion 

The calculated line pr ofiles show some s imilarity with the aver- 
age profile derived in lHonv et al.l (l2009h which peaks at 10 //m. 
The peak position does not match exactly but probably within 
the limit of errors with which such averaged profiles can be de- 
termined. The broad extension to the long wavelength side of 
the observed profile, however, does not agree with the shape 
of our profiles. The strong asymmetry of the observed profile 
is conspicuous. This may have a number of reasons. Either the 
observed profile is the superposition of the strong lOjum emis- 
sion feature of solid silicon monoxide with one or more absorp- 
tion bands of additional dust materials centred around somewhat 
longer wavelengths, or the properties of solid silicon monoxide 
deposited under conditions in circumstellar shells differ from 
that of the material obtained in our laboratory experiments by 
CVD. 

With respect to a contribution by some other materials the 
most likely candidates are hibonite (CaAlgOig) with a strong 
absorption band centred on 12.3 jum, spinel (MgAl204) with a 
strong absorption band centred on 13//m, or corundum AI2O3 
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with a strong and broad band peaking at 1 1 /im. Such mate- 
rials probably would start to condense at higher temperature 
than solid SiO because of their lower vapour pressure. They 
would use up only small amounts of the oxygen because of the 
much lower element abundances of Al and Ca compared to that 
of Si. This would not noticeably compete with SiO condensa- 
tion under conditions where oxygen is scarce. Explorative cal- 
culations showed, ho wever, that the hibonite (optical dat a from 
iMutschke et al.ll2002l) and spinel (optical data from Palikl ll985h 
band would always be seen as a separate band, which seems not 
to be found in spectra of S-stars with a strong lOjjm f eature. The 
rathe r broad corundum feature (optical data from Koi ke et all 
11995) also cannot explain the strong asymmetry of the observed 
profile, but it would fill up the dip of the emission profile due to 
solid SiO at about 12.5 yum. 

Also forsterite may contribute somewhat to the observed 
profile and increase its width and asymmetry, because if the 
oxygen abundance is not close to the critical abundance limit 
eb.crit = fc - fsi + es where O is completely bound by CO and 
by that fraction of Si that is bound in SiO and not in SiS, then 
some oxygen is left over to form silicates. However, it seems 
not to be possible to reproduce the observation with a mixture 
of solid SiO and silicates, because if more than a few per cent of 
the Si forms a silicate, the 1 8 yum feature becomes clearly visible. 
On the other hand, additional formation of silicates besides solid 
SiO could explain the sources with a marked IQfim feature but 
an only weak 18//m feature, which are frequently found within 
the population of S-stars. 

Possibly one could reproduce the observed band structure 
of the stars under consideration by a mixture of the oxides and 
silicates with solid SiO, but this is beyond the scope of this paper. 

Another possibility for the strong asymmetry of the observed 
profile may be that the material formed in stellar outflows is 
more strongly disordered than the material formed in our labora- 
tory experiments. A significantly stronger damping constant yj 
would also result in an asymmetric and much broader line pro- 
file. However, test calculations with strongly increased values of 
jj showed that the observed band shape cannot be explained in 
this way. 

Finally some remarks are in place on the structure of SiO. At 
low temperatures this material slowly d ecomposes into silico n 
nanoclusters embedded in a Si02 matrix (iKamitsuii et alJl2b04l) . 
The material called silicon monoxide then in truth is an inti- 
mate mixture of two different phases which is inhomogeneous 
on scales of a few nanometres. The absorption peak then shifts to 
the characteristic peak of SiOa at 9.2 yum. Clearly, such a decom- 
position did not occur for our silicon monoxide films and would 
also be of minor importance in circumstellar shells if solid SiO 
would be the carrier of the lOyum feature. Decomposition times 
have been measure d for the closely r elated material SiiOs in the 
laboratory by, e.g.. lNuth III & DonnI (fl984 ). If this can be taken 
as representative for solid SiO the decomposition would be suf- 
ficiently slow at the temperatures of interest (< 750 K) that solid 
SiO can survive for a sufficiently long time to be observable. 

5. Conclusions 

This paper reports on laboratory studies of the infrared optical 
properties of solid silicon monoxide. The dielectric function s is 
derived from transmission measurements of thin films obtained 
by evaporation of commercially available solid SiO and deposit- 
ing its vapour on a cold substrate. A Brendel oscillator model is 
fitted to the results of the transmission measurements to deter- 
mine e(w). 



Solid silicon monoxide is a material that has been speculated 
at different occasions to be involved in the formation process of 
silicate dust in circumstellar environments like accretion disks 
around protostars or dust shells around AGB stars. In this paper 
we study the possibility that solid silicon monoxide forms in cir- 
cumstellar dust shells a separate dust component if insufficient 
oxygen is available to build the Si04-tetrahedrons of normal sil- 
icate dust material. This would apply for instance to S-stars with 
a C/O abundance ratio very close to unity. 

The characteristic property of such SiO dust would be a 
broad and structure-less emission feature centred on about 10 /im 
resulting from Si-O bond stretching vibrations, similar to the 
9.7 fim feature of amorphous silicates, but a missing 18// feature 
resulting from O-Si-O bending modes in a Si04-tetrahedron. 
Such peculiar cases have, indeed, been detected for some S-stars 
(Honv et al. 2009). We propose that we see in these stars solid 
silicon monoxide as an abundant dust component. Radiative 
transfer calculations for circumstellar dust shells have been per- 
formed using our new data on the dielectric function of solid SiO 
to calculate the emission band structure due to SiO dust. 

The resulting emission feature at 10//m is co mpared with 
one of the two average emission profiles derived in iHonv et al.l 
(I2OO9I). The profile of the emission feature obtained in the model 
calculation peaks at about just the same wavelength as the ob- 
served feature, which makes its identification as due to solid 
SiO very likely. However, the observed profile is much more ex- 
tended to longer wavelengths than the synthetic profile which 
either may result from blending the band from solid SiO with 
emission bands from a number of additional minor dust species 
(corundum, hibonite, ...) or shows that the laboratory produced 
amorphous SiO films have a somewhat different lattice structure 
than the material condensed in stellar outflows. 

More work is required to explain the asymmetry of the ob- 
served peculiar 10/im feature in some S-stars before the identifi- 
cation of solid silicon monoxide as carrier material of this feature 
can be considered as safe; but in any case this explanation can 
presently be considered as the most likely one. 
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'Deutsche Forschungsgemeinschaft (DFG)'. 



Appendix A: Radiative transfer model 

Model spectra of S-stars are calculated using a simple code for 
modelling radiative transfer in circumstellar dust shells. The ba- 
sic assumptions on the structure of the dust shell follow the con- 
ventional assumptions: A spherically symmetric mass distribu- 
tion around the central star, with inner radius Ri and outer radius 
/?a. The outflow is assumed to be stationary and the outflow ve- 
locity Vexp to be independent of distance r from the centre. The 
radial distribution of mass density in this case is 

M 

eW^-^ . (A.l) 

Here M is the mass-loss rate that is also assumed to be constant. 

With given opacities (see Sect. 14.1b . the radiative transfer 
equation in spherical symmetry 

(A.2) 

or r Oil 
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is solved by the so called p-z- method (cf. lMihaiaslll978h . The 
source function S « and the total extinction coefficient ;h;™' are 



=Z |^^'>(^') + S^ [(3 -^^)A+(3/.^ - 1)K,]^ (A.3) 



The sum runs over all dust species. The angular distribution of 
the scattering term is that for small particles. The temperatures 
r, of the dust species are determined such as to satisfy radiative 
equilibrium 



Jo 



dv;^*^ [Jy - By(Ti)] = 



(A.5) 



for each species by applying a U nsold-Lucy temperature correc- 
tion procedure (see iLucvIl 19 64') adapted to the spherically sym- 
metric case. This requires an iteration procedure for the deter- 
mination of all temperatures Tj. This iteration is combined with 
a simple iteration scheme (successive over-relaxation) with re- 
spect to Jy and Ky in the scattering contribution of the source 
function. 

The inner radius Ri of the dust shell is fixed by the require- 
ment, that the most stable of the dust species appears at some 
prescribed temperature T^. Since it is not known in advance at 
which radius this condition is satisfied, an additional iteration 
procedure is required to determine Ri. The outer radius R.^ is al- 
ways taken at 10^ stellar radii. 

The radiation field of the central star is approximated by a 
black body radiation field. 

Appendix B: Relation to Voigt function 

We show here how the Brendel oscillator model is related to the 
more familiar Voigt function for the profile of damped spectral 
lines. 

A single oscillator of the Brendel model contributes to the 
dielectric function a term 



oo 

TTO-j /' 



IttcTi 



(B.l) 



or, after splitting this into its real and imaginary part, the terms 

0? ,(z2 - a?) 





CO 

Incrj j ^ 



(Z2 - w2)2 _ y2^2 ■ 



(B.2) 



(B.3) 



If 7 <K w, as in our case, the Lorentz profile as function of z is 
different from zero in an only small interval around a>. Letting 
Az - z - CO, \Az\ CO one obtains by retaining only terms of 
lowest order in Az 



CO 

IncTj J ^ 



12a] 



2 '^p,; (Z - CO) 



2w {z - w)2 - iy2 



Ina-j J"'^^ 



4 (z-co)^-iyj 



(B.4) 



(B.5) 



If we compare the imaginary part with definition 



H(a. 



oo 



dxe 



1 



(B.6) 



(x - v)2 + 

of the Voigt function (see ' MihalasI Il978l for details) we rec- 
ognize that the imaginary part of one oscillator in the Brendel 
model corresponds, except for a constant factor, to the Voigt 
function with parameters 



a - 



CO - coo 
V2cr 

8cr ■ 



(B.7) 
(B.8) 



Since the Voigt function is the real part of the complex er- 
ror function (e.g. Humlicek 1982, who also gives a FORTRAN 
routine for calculating the complex probability function), the 
Brendel oscillators are also directly related to the complex prob- 
ability function. 
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